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Abstract 

Incorporating all recent theoretical advances, we resum soft-gluon corrections to the total tt cross-section 
at hadron colliders at the next-to-ncxt-to-leading logarithmic (NNLL) order. We perform the resummation 
in the well established framework of Mellin A^-space resummation. We exhaustively study the sources of 
systematic uncertainty like renormalization and factorization scale variation, power suppressed effects and 
missing two- and higher-loop corrections. The inclusion of soft-gluon resummation at NNLL brings only 
a minor decrease in the perturbative uncertainty with respect to the NLL approximation, and a small 
shift in the central value, consistent with the quoted uncertainties. These numerical predictions agree with 
the currently available measurements from the Tevatron and LHC and have uncertainty of similar size. 
We conclude that significant improvements in the tt cross-sections can potentially be expected only upon 
inclusion of the complete NNLO corrections. 



1. Introduction 



The production of tt pairs at hadron colliders is well understood within next-to-leading order (NLO) 
perturbative QCD, where corrections of order ©(af) are included. Results have been available for a while 
for the fully inclusive [l], 0, one-particle inclusive B lU, two-particle inclusive production H, including 



decay Ml, spin correlations [8|, |9j], off-shell effects [lO|, and associated jet production [12, ll3|, llJ, ll5 
Til . ItI [Tg ■ In addition, the resummation of next-to- leading logarithmic (NLL) soft gluon effects has been 
lon g es tabl ished Il9i. and results beyond the NLO and NLL level of accuracy have recently been published 



[201 l2l|, l22|, [23], [2 ^ 

At the current level of precision, the theor y ag rees 20l 23 . 26l 28 . 2^ with the data from the Tevatron 
and LHC [13, [3l|, [11 111, [11 [11, [11 113, S 111, Efl The large statistics becoming available at the LHC, 
and the precision of the experiments, will however soon bring the accuracy of the measurements to the level 
of possibly less than 5 percent, thus challenging the present theoretical systematics. 

In this paper we extend to the ncxt-to-next-to-leading logarithmic accuracy (NNLL) the Mellin A^- 
space soft-gluon resummation of the top-pair production cross-section which was first developed in (4^ |43| 
for leading (LL), and in [l^ for next-to- leading logarithmic accuracy (NLL). This calculation uses the 
recently derived two- loop anomalous dimension matrices [H EBj , and matches the result to the best known 
approximation to the NNLO cross-section [22j- The numbers we derive are very robust and show a significant 
stability - as expected - given that this is the third order at which the soft-gluon resummation is applied to 
this observable (i.e. LL, NLL and NNLL). 

We include in our assessment the latest sets of parton distribution function (PDF) fits, analyze the 
impact of the emerging contributions at 0{a1) like two-loop Coulomb corrections and, most importantly. 
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missing NNLO corrections. Addressing the full theoretical uncertainty associated with the total tt cross- 
section as currently known, we conclude that we see no evidence for a strong reduction of the theoretical 
uncertainty compared to the long-ago established NLO-I-NLL analysis. Based on our comprehensive analysis 
we, however, speculate on the uncertainty that should be achievable once the full NNLO calculation will be 
completed, which could be significant. 

The paper is organized as follows: in section [5] we summarize the key elements of the threshold approx- 
imation at NNLO, and develop the Mellin-space resummation formalism to the NNLL level, leaving some 
technical details to the Appendix. In section [3] we study the details of the theoretical systematics and the 
sensitivity of our theory prediction to the scale dependence and to a number of currently unknown contribu- 
tions of 0{a*) and higher. In section |4] we include the study of the PDF systematics, and present our best 
predictions for the Tevatron and LHC (7 TeV). We present the dependence of the predicted cross-section on 
the value of the top mass and compare with the most precise available experimental measurements, finding 
a very satisfactory agreement. In the concluding section we discuss the comparison of our results with the 
current literature. 



2. Theoretical framework 



2.1. NNLO results in the threshold approximation 

While the calculation of the full NNLO results for the total cross section is yet to be completed, recent 
work has determined its exact behavior near the production threshold, where it can be represented as an 
expansion in powers of 1//3 and log (3 {f3 ~ ^^1 — 4m^/.s being the heavy quark velocity in the QQ rest frame 
and m the top pole mass). In the limit /3 — >■ the 0{a'^) partonic cross section can be written as 
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where we defined al^\ as the Born- level cross section. For ease of notation we set the factorization {^ip) and 
renormalization {(ir) scales equal to /i, although in our calculation we allowed them to vary independently. 

The index 1=1,8 corresponds to the color configuration of the heavy quark pair, (ij) ~ {qq), (gg) are 
the initial-state partons and the functions fjf^.n) ^^.^ expanded as follows: 
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For n ~ 1,2 the above functions, including the /3-independent terms C^^'"\ can be derived entirely from the 
knowledge of the NLO results. In particular, the C*-^'"^ (n ~ 1,2) terms can be determined by enforcing 
the compensation of their scale dependence against that of appropriate terms of ©(a^). For n = 0, the 
functions were extracted in from the threshold behavior of the NNLO result, up to (but excluding) the 
constant term C'^'''-', which will only become available with the completion of a full NNLO calculation. 

To the extent that the hadronic production of heavy quark pairs receives an important contribution 
from the region near threshold, it is meaningful to incorporate the known singular threshold terms into an 
improved prediction of the total production cross section, even though they do of course not represent the 
full NNLO resuh. We label NNLO/j [H] the approximation that adds to the fuU NLO ah 0(af ) terms that 
become singular when /3 — J- 0. Note that we do not impose scale dependence cancellation at ©(a^), i.e. 
we exclude all the finite contributions C^^'"': while the n = 1,2 terms are known from scale dependence 
compensation with ©(af) terms, the lack of any threshold enhancement prevents us from assuming that 
their value should be bigger than the unknown C'^'"-' (This will be confirmed with direct numerical studies 
presented in Section [3.21 ). Thus we find it more coherent to neglect the non-singular terms. We also 
remark that if one wants to include these scale dependent terms, one should still consider the ambiguity 
in the choice of the scale that divides the renormalization and factorization scale in the arguments of the 
logarithms. These scales are uniquely fixed only if the constant terms are known. Thus, the compensation of 
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the rcnormalization and factorization scale variations induced by these terms wiU be counterbalanced by the 
uncertainties obtained by varying these new scales. It is therefore simpler not to include these terms at all, 
and let the lack of scale compensation work as an indicator of unknown higher order terms. Our approach 
is therefore to account for the ignorance of the 0{a'^) constants through the scale variation systematics. In 
addition, we shall show explicitly in the following that two different choices of these constants, C**^^'"^ = 

and C*^^'"-* = C"^^'"' (see Appendix), lead to differences consistent with the scale systematics. 

2.2. Soft-gluon resummation with NNLL accuracy 

We extend here the NLO+NLL formalism for soft gluon resummation of the total top-pair hadroproduc- 
tion cross-section, discussed in Ref. [l^, to include the resummation of NNLL terms. 
For simplicity, we first summarize the qualitative overall structure of our result: 

• the 0{a\) contributions correspond to the exact NLO result; 

• We perform the resummation in Mellin space and then invert numerically back to a;-space using the 
Minimal Prescription of Ref. [4^ . 



the truncation of the NNLL result at 0{oii) includes all singular contributions described by NNLO^, 
plus non singular terms that arise from the inverse Mellin transform of the A^-space resummation; in 
particular, terms of 0{\/N) can give non-negligible contributions, which reflect the uncertainty about 
higher-order non-singular terms. This point in particular underscores a qualitative and quantitative 
difference with the alternative approach of simply matching at 0{al) the resummed result to NNLO^. 

The resummed cross-section in A-spac^B reads: 



[res) 



v=<]<],gg 



+ -T''\r,^^) , (3) 

where cr^'^'^'' is the NLO hadronic cross-section ^Ij'^n NNLL resummed partonic cross section, 

and (ctI^%''i is its perturbative truncation at order a^. As before, we set here the rcnormalization (/i/?) 

and factorization (hf) scales equal to fi: they are however kept separate and varied independently in the 
subsequent studies of the scale systematics. 

In the threshold limit N ^ oo the NNLL resummed partonic cross section ct^J']^'' factorizes: 

/ 2 

-(res)/ 2 r 2\ 2\ ^(Coul), , 2\ 2 / 2\ -{Hard), / 2\ 2 i 2\ j\ ( / 2\ 

<^lj,Ni'^ '"4/^ ) = 2^ %m,a/("4M ),fJ- /m ) [asifJ- ),fi /m ) Ay,i,7v+i I ), ^ 
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in terms of the radiative factors n containing all contributions due to soft-gluon emission, the func- 
tions o'l^j"^'' containing the threshold-enhanced bound-state contributions (Coulomb effects) and the hard 

matching functions a\^i^'^\ 

The radiative factors introduced in Eq. A,;jj_Ar, are given by: 

In Ay- 1, jv = / dz <^ / 25,,A, {as{ql)) + Aj.i {a, (477x^(1 - zf)) \ 

+ 0{as{as\nNf) . (5) 



-"^The Mellin moments of a function g{p), with p = Am? /s, are defined by gjy = dp ^ g{p)- 
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The process independent anomalous dimensions Ai describe soft-coUinear initial state radiation. They 



are known through three loops [47|, |48|, |49[ ; the explicit expressions can be found in Ref . [49| . The anomalous 
dimensions Dij i describe wide-angle soft radiation and depend both on the initial and final states [H^ . For 
top- pair production, they have recently been derived through two loops in Refs. [H, E^- Their explicit 
expressions read: 






108 36 ' y ^ V54 

where Np is the number of light flavours, i.e. lighter than the top quark. The integrations in Eq. ([S]) can 
be performed analytically and the resummed NNLL cross-section can be written explicitly. To extend the 
NLL results of Ref. [l^, one can utilize, for example, the results in Ref. [5l| . 
The A'^-independent, hard matching functions read: 

(7) 

The one-loop results H^j^\{ii^ / m?) have been calculated in Ref. [H^l- The renormalization/factorization 

scale dependence of the two-loop corrections h\J'^{^'^ / m?) can be obtained, for example, from the results of 
Ref. [13: 

^^SsCmV™') = ^^8(1) + 8.91918 In^ (^^^ + 34.7212 In (^^^ , 
^SsCm'M') = ^Ssll) + 9.31619 In^ + 60.2080 In 



Hgf^iifJ'Vrn'^) - ^Si(l) + 9-31619 In^ + 38.5239 In ) . (8) 

f 2 i (20) 

The genuine two-loop constants -ffj-j i(l) arc currently unknown, and are related to the constants C*,-^ j 
introduced in 1^, as discussed in the Appendix. For simplicity, in Eq. ([5|) we have given directly the 
numerical values of the scale dependent term evaluated for Np — 5. In analogy with our discussion of the 
finite C*^^'"-* coefficients, which we suppress in the NNLO^ approximation, we shall set ij|Jj(/i^/TO^) = in 
our default NLO-I-NNLL predictions, as will be motivated by the numerical results of Section [321 

To analyze the impact of subleading, power-suppressed corrections we follow Ref. [l^ and introduce 
a power suppressed term controlled by a constant A into the function H^l\. The choice A = sets this 
additional term to zero, while the default value A = 2 was chosen in Ref. [19[. The rationale behind the 
inclusion of this term was the observation that power suppressed terms 0{1/N) are needed to bring the 
resummed cross-section closer to the fixed order NLO result away from the threshold region. We will have 
more to say about the numerical impact of this term in Section [31 

In deriving the Coulomb contributions o'j^j"^'' we follow the approach of Ref. [19} and absorb in it the 
exact Born cross-section, except for its overall factor of that we attribute to the hard function - see Eq. ([7]). 
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The functions o'^^^"!^ {as{^\), ^j^/m^), which only depend on the renormahzation scale, are obtained as the 
Mellin transform of the following x-space functions: 
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The Born cross-sections, for all reactions and color configurations, have been given in Ref. [19[, together 
with the one-loop Coulomb functions C^^^if) = Cpn^ /{2P)and C^j^ip) = {Cp - CA/2)7rV(2^). The 
two-loop functions can be extracted from the results of Ref. [2^ by matching them to the Mellin-inverse of 
Eq. 
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The Mellin transform of Eq. © through NLO has been calculated in Ref. [19[. We do not present here 
the Mellin transform of the NNLO corrections. They are given by large expressions that are straightforward 
to calculate following the discussion of Ref. [l9l | and utilizing the following approximation: 
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Before concluding this section we stress again that in Eq. ([3]) we use the NLO fixed order result 
^(m^o)(^^'^ and not cr(^^'^'-''^)(TO^) as one might expect. The reason is that all the information to be 
found in the approximate NNLO cross-section cr'-'^^'"'-''*^ is already contained, hy matching, in the all order 
resummed result o'f^^ . In particular, the two loop anomalous dimensions in Eq. ([5]) control the single In /3 
terms in cr^'^^^°''\ and the two-loop Coulomb terms in Eq. ([3]), including the potentials ~ ln/3, have been 
matched to the Coulomb terms in (t'-'"'^'"'-'^'^. 

To make this point completely transparent we note that the difference between Eq. ([3]) and its analogue 
defined by using cr'-^^'"'^'^' instead (and, of course, subtracting the terms in of-^^ through 0{oi\)) is given 
in N-space by the following expression: 
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In other words, Eq. ([T^ represents the difference between the terms of order derived respectively within 
the resummed A^-space and the fixed order x-space approaches. Since, as we just explained, the two contain 
the same input they do cancel each other, at least in the limit N 00. In practice Eq. contains 
power-suppressed terms that behave as 0{1/N) in the soft limit iV — > 00. These power suppressed terms 
originate in the lower loop (LO and NLO) terms in Eq. (j4] ) and, as it turns out, are not numerically 
negligible. Given that both terms in Eq. (|12p are the result of an approximation, and as we favour TV-space 
resummation since it's less likely to introduce large terms due to the violation of momentum conservation, 
we perfer not to introduce the power terms Eq. into Eq. ([3]). Such 0{1/N) ambiguity is inherent in the 
soft approximation irrespective of the details of its implementation and can only be removed by adding to 
Eq. ([3]) the full NNLO result, once it becomes available. 
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3. Study of the theoretical systematics 



In this section we focus on the purely theoretical systematics, arising from the scale dependence of 
the cross sections, and from the different possible descriptions of higher-order terms not controlled by 
resummation, like unknown two- loop (threshold) hard matching constants ^jj i(l) and terms vanishing 
at threshold. We shall then complete the study of systematics, including PDF and mass dependence, in 
Section m 

3.1. Benchmark results 

We shall compute reference values for our ti cross section predictions at mt — 173.3 GeV (STipl. The 
central values of these predictions are obtained for /i^j = fip ^ mt. Throughout the paper we use the strong 
coupling constant evaluated at scale fj,ii as provided by the corresponding PDF set. Our default parton 
distribution set for NLO (and NLO-^NLL) is MSTW2008nlo68cl, whereas for NNL0;3 and NNLL resummed 
calculations we use the MSTW2008nnlo68cl set [H^. In all cases we include them through the LHAPDF 
interface (HTj . 

The scale systematics is evaluated by varying the renormalization and factorization scales independently 
in the range suggested in ref. fisj : 

mt/2 < fiR, fip < 2mt, with 1/2 < fiji/ fip < 2 , (13) 

and searching for the minimum/maximum of the resulting cross-section. It is usually sufficient to consider 
only the endpoints of the range Eq. ([T^ . namely the pairs {^r/mt, np/mt) = (2,1), (0.5,1), (1,2), (1,0.5), 
(2, 2) and (0.5, 0.5) . We have verified that a search over a grid with a few hundred points satisfying Eq. (|13p 
in the (^FjA'a) plane agrees to within few per mille with the minimum and maximum rates found in the 
scan of the endpoints. 

Different power-suppressed terms are probed by varying the parameter A over the two values ^ = and 
A — 2, as discussed in [19[. 
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Table 1: Central values and theoretical systematics for the various approximations to crtot, in pb, at the Tevatron. mtop = 173.3 
GeV, PDFNLO=MSTW2008nlo68cl, PDFNNLO=MSTW2008nnlo68cl. Row 9, highlighted in blue, gives our best prediction 
for central value and scale systematics. The predicted cross-sections are presented, if applicable, depending on the values of 
the constant A and on whether pure two-loop Coulomb corrections JTUJ are included or not. 



^The best measured value for the top mass has been recently updated in Ref. [53^ to 173.2 it 0.9 GeV (for a recent review 
see [541). However, for use here the previously published value of 173.3 to facilitate the comparison with other recent theoretical 
analyses that used this nit, such as [55^ . We estimate (see section |4]l that the change of 0.1 GeV, from 173.3 to 173.2 GeV, 
would lead to an increase of about 0.3% in the cross section values, well within the overall theoretical uncertainties. 
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Table 2: As in Table[T]but for the LHC at 7 TeV. 



Our numerical results are summarized in tables [T] and O As a benchmark, in rows 1-3 of these tables 
we present the well-understood NLO and NLO+NLL {A = and A = 2) results. They are an update of 
Ref. 0. 

In row 4 we present the NNLO^ approximation as defined in Section 2, and in row 5 we show the effect of 

including non-zero values for the constants C^J'"'', setting them to the value c[^'^^ defined in the Appendix 
in Eq. (jA.2|) . We notice that, while the numerical impact of including these constants is noticeable, it is 
smaller than the overall scale uncertainty. Therefore we argue that the scale variation can largely account 
for the uncertainty stemming from unknown part of the higher-order terms. 

At order the fixed order NNLO approximation contains terms that are not predicted by the NLO+NLL 
result, like NNLL soft-enhanced terms and pure two-loop Coulomb terms, i.e. Coulomb terms that do not 
arise in the expansion of Eq. Q from the product of one- loop contributions to aj^""'-' with ct^^"'"'''' Ajv+i- 
It lacks however terms at O(a^) and beyond that are contained in the resummed results. The NLO+NNLL 
approximation combines both these ingredients, and is superior to the NNLO^ one, since it contains all the 
information to be found in NNLO^, plus the towers of soft LL, NLL and NNLL logs beyond order 0(q!^). 
The NLO+NNLL rates are given in rows 6-9, where we also describe the impact of subleading terms 
(through the constant A), and of the two- loop Coulomb effects that were absent in the NLO+NLL results 
of Ref. 0. 

We observe that these two-loop Coulomb terms ([TO)) have a sub-per-mille effect on the central values 
for the Tevatron. At the LHC the effect is larger, of order 2%. Their effect on the scale uncertainty is at 
most at the few per-mille level and is thus negligible (for both Tevatron and LHC). We do not resum the 
Coulombic corrections beyond order aj.. This resummation has been performed in Ref. [HBI and the effect 
was found to be negligible. 

From tables [T] and [5] we also observe a dependence of the predicted cross-section on the value of the 
constant A. As we emphasized before, we consider the inclusion of non-zero A as a model for the power 
suppressed terms ^ that are not controlled by the threshold approximation. We observe a modest 0.5% 
decrease in scale dependence from including A and a 2% shift of the central value, which is consistent 
with the overall scale systematics (the corresponding changes for the LHC are 0.5% and 3% respectively). 
We conclude that, while our estimate of the size of the power suppressed terms is not comprehensive, it 
clearly shows that power suppressed terms can be a significant, few-percent-effect on the central value both 
at the Tevatron and LHC. 

Before closing this section we offer an alternative graphical representation of the scale uncertainty of 
the various scenarios given in tables [T] and [51 In Fig. [T] wc plot the relative half uncertainty due to scale 
variations (defined as Adscaics/c = (cmax — crmin)/(2crcontrai), so that one can quote for the cross section 
^central ± AcTgcaics) 8-8 a function of the hadronic centre of mass energy for pp collisions and for a top mass 
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Figure 1: Relative half uncertainty in the total inclusive tt cross-section as a function of the centre-of-mass energy V-S'had of 
a proton-proton collider. Left plot: whole energy range up to V'S'had = 50 TeV. Right plot: blow-up of the threshold region, 
up to y/S^ = 2 TeV. 



of 173.3 GeV. As expected, the uncertainty of the NLO calculation is largest when the production is closest 
to threshold, and decreases for larger -v/S'had- Upon resummation of the threshold logarithms, the biggest 
improvement (i.e. reduction in uncertainty) can likewise be obtained close to the threshold. One can easily 
observe the better stability of the NLO-I-NLL and NLO-I-NNLL results when the threshold is approached. 
As the ti production takes place at larger ySh^ the effect of resummation is reduced, as expected, since the 
rcsummed logarithms become smaller. When Shad 3> 4m^ one does not expect any significant improvement 
from a resummation of logarithms that are strongly suppressed. Indeed, one observes from the plot that 
the uncertainty of the rcsummed results is practically identical to that of the fixed order calculation in this 
limiiH. 

3.2. Impact of scale-dependent finite terms at ©(af ) 

We discuss here the effect of the 0{a'^) constant terms C^-^'^\ which wc introduced in ^ and which 
we suggested should not be incorporated in either NNLO^ or in the NNLL rcsummed results. Their main 
impact is the large reduction in the scale variation of the cross section, due to the scale logarithms present 
for n = 1, 2, whose coefficients compensate by construction the scale dependence of ©(af) terms. If we focus 
on the NNLO^ case, the addition of the known ClJ'"' contributions {n = 1,2), for different choices of the 
unknown n = terms, leads to the following results: 

(14) 
(15) 
(16) 
(17) 

The central values coincide with those obtained in absence of the clj'""* (n — 1, 2) terms, since the logarithms 
vanish at the central value /x = m. However, the scale dependence is much smaller than in rows 4 and 5 of our 
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, r. f, + 2.8(1.8%) , 
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for 


/^(2.0) _ 






^The NNLO^ result seems to display a slightly smaller uncertainty than the resummed ones. The difference is however 
likely not significant and, in particular, it should not be considered as suggestive that this approximation constitutes a better 
prediction. 
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previous tables. At the Tevatron (LHC), the scale dependence drops from about ±5% (±7%) to about ±3% 
(±2%). This significant reduction clashes however with the comparable or larger cross section variations 
(3% at the Tevatron and 4% at the LHC) induced by the variation of cfj'°'' within an a priori reasonable 

range. We conclude that the significant reduction in scale dependence in presence of C|J'°'' (n = 1, 2) terms 
cannot be interpreted as a genuine reduction in the total theoretical uncertainty, unless it is combined with 
the systematics emerging from the unknown value of the two-loop constants C^j'^K 

Similar conclusions can be drawn from the study of the NLO+NNLL results. From the viewpoint of 
threshold expansion, in the limit — > oo one has to neglect the unknown constants H^j'^^l) and the A^- 

independent, fxp, //ij-dependent logarithmic terms, i.e. the whole two- loop hard function H^'^\{fi'^ / rrt^) . This 
is completely analogous to what happens in the NNLO^ approximation. In case the scale dependent terms in 
H^'^\{^'^ /m?) are retained, the theoretical uncertainty should include the variation of the unknown two-loop 
constant again as in the NNLO^ approximation. For the sake of documentation, we quote here 

the relevant results for the scale dependence obtained after inclusion of the hard function ij|Jj(//^/m^), 
exploring as an example the two cases of hIJI{1) = and -ff^J^i(l) = H^-^lil) introduced in the Appendix: 

Tevatron: a^LO+NNLL"'^' ^ g 722+0;Oi^™ pb ^(J) (i) = q (18) 

^NLO+NNLL«-^g,9g3+0.009(a^^^ H^^il) = 11^(1) (19) 

LHC : ^NLO+NNLL«'=' ^ ^53^+5.6(3.6%) j^^. ^(2)^^^ ^ q ^20) 

^NLO+NNLL«<^'^,g7 9+5.2 (3.1gp^ for HglH) ^H^,(l) , (21) 

where A = 2 throughout. These scale uncertainties are slightly larger than those found for NNLO/3 in 

_____ 

([T ^ - pT]) . but they are still small compared to the impact of the unknown finite contributions of i?jj i(l), 

as suggested by the comparison between the H ~ Q and H = H results above. 

(2) 

These remarks justify our choice not to include the H^^ j function in our benchmark predictions for 
the central value and the theoretical systematics. We notice nevertheless that the significant reduction in 
scale variation obtained with the inclusion of the known, finite scale dependent terms, is indicative of the 
uncertainty of the full NNLO result, when it will become available. 



4. Phenomenology 

It is well known that the purely theoretical uncertainty related to the lack of high-order corrections is 
only a fraction of the overall systematics. Recent studies of the tt cross section uncertainty due to the PDF 
parameterization, including the latest fits by several groups, have been reported in |5 a. | 59| . These studies 
considered the fixed-order NLO results and the approximate NNLO calculation of Ref. [20[, as implemented 
in the program HATHOR [6^. The main conclusion of those studies was the consistency, at the LHC 
energy, between the central values and the uncertainty bands obtained using most PDF fits (MSTW08 [H^ , 
NNPDF2.1 [59II, GJR [6l|, CTIO both at NLO and NNLO, with some differences, incompatible with 

the estimates of the systematics, with respect to other sets such as ABKM09 [6^ and HERAPDF [63 |. 

Given the minor changes in fixed-order versus resummed predictions, we expect that the PDF uncer- 
tainties estimated in 5^ 5^ should not be affected by resummation. We verify this result explicitly here, 
considering our default PDF set MSTW2008nnlo68cl and the PDF set NNPDF21_nnlo_nf5_100 [Hij. 

Our results for the total tt cross-section using MSTW2008nnlo68cl are: 



NLO+NNLL 
'tot 



(Tevatron; rrit = 173.3 GeV) 



6.722 +°-!fJ^lS [scales] tlf.Z^^l [PDF] pb , 



-0.410(6.1%) 



NLO+NNLL 
"tot 



(LHCttov; mt = 173.3 GeV) 



158.7 



+ 12.2 (7.7%) 
-13.5(8.5%) 



[scales] 



+4.3(2.7%) 
-4.4(2.8%) 



[PDF] pb , (22) 
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Figure 2: Theoretical prediction for the total inclusive tt cross-section at the Tcvatron, as a function of the top mass, versus the 
measurements of Refs. [S^l , [H . The left plot is for the MSTW2008nnlo68cl PDF set, the right plot for NNPDF21_nnlo_nf5.100. 
The uncertainty is a linear sum of scale uncertainty (the white central band) and PDF uncertainty (red bands). The central 
value is shown with a black line. The theoretical predictions in this figure correspond to row 9 in tabled The horizontal bars 
on the measurements reflect the uncertainty in the measured top mass (see footnote [2] ) . 



where we defined the upper and lower limit of the scale variation using the endpoint scan defined after 
Eq. (O. With NNPDF21_nnlo_nf5_100 we obtain instead: Q 

^NLO+NNLL(Tevatron; = 173.3 GeV) = 7.021 ttllllli^i ^es] tlf^Z'^^J, [PDF] pb , 

^NLO+NNLL(LHC,Tcv; m = 173.3 GcV) = 163.1 tllZ'vi ttZ-l^], [PDF] pb. (23) 

The PDF uncertainties in both these sets of predictions should be considered to be at the 1-a level. 

A few comments are in order. To start with, we confirm that the relative PDF uncertainty at NLO+NNLL 
is similar to that derived from a fixed order calculation, as in [H^, [H^ . We also notice that the scale uncer- 
tainty is rather independent of the PDF set. This is consistent with the fact that the relative contribution 
oi gg, qq and qg initial states does not change significantly when changing PDFs. 

We also confirm the consistency of the central value and of the PDF uncertainty estimated, for the LHC 
at 7 TeV, using the MSTW and the default NNPDF2.1 sets. We note, on the other hand, that at the 
Tevatron the NNPDF prediction is larger than MSTW by about 5%, compared to the individual estimates 
of PDF systematics, which are of the order of ±2%. This difference can be understood in terms of the 
different values of the strong coupling constant associated with the MSTW2008nnlo68cl {as{Mz) = 0.117) 



and the NNPDF21_nnlo_nf5_100 {as{Mz) = 0.119) sets, and the fact that the cross section scales like a^. 
To better quantify this effect we have also used the set NNPDF21_nnlo_as_0117_100 [H^ to compute the 
central values corresponding to Eq. (^5)) : 

fTf„^°+^^^^ (Tevatron; mt = 173.3 GeV) = 6.742 pb, 

^nlo-hnnll^LHCttcv; "It = 173.3 GcV) = 156.8 pb. (24) 

These results confirm that the apparent inconsistency between the Tevatron predictions of the NNPDF and 
MSTW NNLO sets disappears when they both use the same coupling constant {as{Mz) = 0.117 in this 
case). This also suggests that an additionally uncertainty of the order of ± 1-2% should likely be added to 



^As a central value we take the number derived with the central NNPDF set, not the mean over the whole set of PDFs. 
The difference is at the per mille level and thus completely negligible. 
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Figure 3: Theoretical prediction for the total inclusive tt cross-section at the LHC @ 7 TeV as a function of the top 
mass versus the measurements of Refs. [l^l.El- The left plot is for the MSTW2008nnlo68cl PDF set, the right plot for 
NNPDF21_nnlo_nf5_100. The uncertainty is a linear sum of scale uncertainty (the white central band) and PDF uncertainty 
(red bands). The central value is shown with a black line. The theoretical predictions in this figure correspond to row 9 in 
table |2] The horizontal bars on the measurements reflect the uncertainty in the measured top mass (see footnote [2] ) . 



any tt total cross section evaluation as a result of the uncertainty with which is known, if not already 
included with the PDF uncertainty. 

We plot in Figs. [2] and [3] our predictions for the total tt cross-section as a function of the top mass in 
the ranee 168 - 178 GeV for both Tevatron and LHC {VS = 7 TeV). In view of the difference between the 
MSTW and NNPDF results for the Tevatron, we present the results for the two PDF sets on different plots. 
The uncertainties from scales and parton distributions quoted in Eq. (j22[) are added linearly. On the same 
figures we compare our theoretical predictions with the most accurate available experimental measurements 
from the Tevatron [s^, [s^ and LHC [4^,[4lj. We display the experimental points at the current world- 
average value of rrit = 173.2 ± 0.9 GeV, without applying any correction factor to account for the difference 
in experimental acceptance with respect to the nit values used in the measurements. These amount to a 
reduction in rate at the sub-percent level. 

We observe that the uncertainties of the current theoretical predictions and experimental measurements 
are comparable in size. The predictions agree with all measurements within the uncertainties, although at 
the Tevatron the data tend to be on high side of the theoretical band, particularly in the case of the MSTW 
cross sections. 

For ease of use, we have fitted the mass dependence of the ^^^0+^'^^^ predictions relative to the MSTW 
sets using the functional form 

aim) = alrriref) { -] 1 + ai ^ + 02 -] • (25) 

\ m J y TUref V "n^ref J J 

The resulting parameters, for the central curve as well as the scales and PDF uncertainties separately, are 
collected in table [S] They provide fits that are accurate to within about one per mille in the 150-200 GeV 
mass range, but should not be used indiscriminately much beyond this region. 



Finally, due to the interest in current experimental searches of possible [65| fourth-generation quarks, both 



at the Tevatron [btI . IgsI I and at the LHC [69l | , we extend our results to the production of a hypothetical 
very heavy quark, T. Its production cross section at the LHC (for ^/S = 7 and 8 TeV) up to my — 1200 
GeV is shown in figure |4l and the detailed breakdown of the systematics is given, for a set of mass values, 
in Table U) The cross sections are calculated by simply changing the mass parameter value in the top cross 
section calculation, thus neglecting the small corrections due to the top quark in the evolution of and of 
PDFs. 

The overall uncertainty (scale-|-PDF) is roughly uniform at the ±10-15% level in the range of the plot. 
This is the result of a decreasing scale uncertainty, which is more than compensated by an increasing PDF 
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iriref = 173 GeV 


a{mref) (pb) 


ai 


02 




Central 


6.785 


-1.394 


7.451 xlO-i 




Scales + 


7.030 


-1.409 


8.047 xlO-i 


Tevatron, pp at — 1.96 TcV 


Scales — 


6.370 


-1.379 


6.919 xlQ-^ 




PDFs + 


6.946 


-1.373 


7.106 xlO-i 




PDFs - 


6.669 


-1.408 


7.527 xlO-i 




Central 


160.1 


-1.191 


8.042 xlO-i 




Scales + 


172.4 


-1.224 


9.096 xlQ-^ 


LHC, ppat ^=7 TcV 


Scales — 


146.5 


-1.162 


7.957 xlO-i 




PDFs + 


164.4 


-1.175 


7.867 xlO-i 




PDFs - 


155.7 


-1.205 


8.416 xlO^i 




Central 


228.6 


-1.069 


6.798 xlO-^ 




Scales + 


246.8 


-1.104 


7.335 xlO-i 


LHC, at Vs = 8 TeV 


Scales — 


208.8 


-1.042 


6.299 xlO-i 




PDFs + 


234.2 


-1.054 


6.533 xlO-i 




PDFs - 


222.7 


-1.083 


6.964 xlQ-i 



Table 3: Parameters resulting from the fit of the functional form in eq. 12511 to our best prediction, NLO-I-NNLL, for the top 
cross section at the Tevatron and the LHC (7 and 8 TeV) as a function of the top mass, including the uncertainties from scale 
variations and PDFs (the MSTW2008nnlo68cl set). These parameters provide fits that are accurate to within about one per 
millc in the 150-200 GeV mass range, but should not be used indiscriminately much beyond this region. 



uncertainty, due to the large partonic x values probed by the production of a large mass object, and PDFs 
being generally less well known in this region. 

Finally, we present our best predictions for the LHC configuration foreseen for the 2012 data taking 
{y/S =8 TcV), with the PDF sets MSTW2008nnlo68cl 0: 

^NLO+NNLL(LHC8Tcv; = 173.3 GeV) = 226.6 tlllfslyj) [^^^^M -5:8(2:6%) PDF] pb . (26) 
and NNPDF21_milo.nf5.100 

^NLO+NNLL(LHC8Tcv; m, = 173.3 GeV) = 233.5 ^^as^sS I^^^M -6:5(2:8%) PDF] pb . (27) 

All of the numerical results presented in this section, and their extrapolation to different values of the 
heavy quark mass or of the scale parameters, can be obtained through a simple web interface [t^. These, 
and more general results obtained under the various approximation scenarios outlined in Section [31 can 
furthermore be computed with the help of the program Top++ fn\ . 



5. Concluding remarks 

Using recent theoretical developments, we extend in this paper the soft-gluon resummation of the total 
tt cross-section to the NNLL order, using the Mellin A^-space formalism. The result includes all known 
NNLO terms that are singular at the production threshold. The current work represents the third-order 
logarithmic improvement for this important collider observable that has been instrumental in developments 
in precision collider physics. 

We explored the implications of the NNLL approximation in a comprehensive phenomenological study of 
the total cross-section, to quantify the full theoretical uncertainty currently associated with this observable. 

In fixed order calculations the theoretical uncertainty is typically identified with the residual scale sensi- 
tivity of an observable. While not perfect, such a procedure is well understood and gives a meaningful way of 
comparing theoretical predictions across different observables and levels of precision. The procedure relies on 
the following considerations: since the exact result must be scale-independent, and the scale dependence of 
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Figure 4: Theoretical prediction for the total inclusive production cross-section at the LHC (7 and 8 TeV) of a TT pair of 
hypothetical heavy fourth-generation quarks. Effects from the running of top quarks in loops are expected to be very small 
and are neglected in these numerical predictions. 



each fully calculated order must be of higher order, one assumes that the residual scale dependence of a cal- 
culation is numerically comparable to the higher-order scale-independent terms, which can only be obtained 
via the complete calculation, and whose size the theoretical systematics attempt to estimatclfl. When dealing 
with an approximate NNLO calculation for tt hadroproduction, and trying to assess its uncertainty via the 
residual scale dependence, one must keep in mind that known terms of 0{a^) include: (a) terms singular 
at the production threshold, both scale dependent and independent, whose behaviour at higher orders is 
determined by general dynamical considerations, and which can therefore be included and resummed, with 
a genuine improvement of the accuracy; (b) finite, but scale-dependent, terms, whose value can be fixed by 
imposing full 0{a'l) scale-independence. Inclusion of such terms will lead, by construction, to a reduction of 
the scale dependence, but this reduction does not reflect the real size of the theoretical uncertainty, which 
is rather governed by the unknown constant terms of ©(af). 

In this paper we assessed the possible size of several unknown higher-order contributions, and studied 
their contribution to the theoretical uncertainty. In particular, we demonstrated that the reduced scale 
sensitivity, obtained by using the exact C(ag) scale dependence, leads, once the uncertainty of the unknown 
terms is accounted for, to a larger overall systematics, comparable to that of the NLO-I-NLL cross section. 

We also demonstrated that the predicted cross section has a few-percent sensitivity to currently unknown 
1 /N suppressed terms that are beyond any of the approximations available in the literature. Summarizing 
these observations we conclude that at present the total uncertainty of the total ti cross-section at the 
NLO-I-NNLL order is only modestly lower compared to the long-established NLO-I-NLL result. Guided 
by the small scale dependence of the results obtained imposing the exact C(q:^) scale dependence, we 
nevertheless speculate a significant decrease of the theoretical uncertainty in the total ti cross-section once 
the full NNLO result becomes available. 

Finally we would like to briefly compare our work with theoretical works that have appeared in the 
recent past and that make, to various extent, use of NNLO approximations. 



^ A detailed presentation of these well-known assumptions is given in 1721 , where they are then used to argue for an alternative 
way of characterizing the pcrturbativc theoretical uncertainty. The method proposed in this paper has however so far only 
been detailed for e+e~ collisions, and cannot therefore be applied to hadronic top production. 
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LHC, ^fS = 7 TcV 


LHC, ^^S = 


8 


TeV 




TUT (GeV) 


o-tot [pb] 


Scale 


PDF 


o-tot [pb] 


Scale 


PDF 


200 


74.71 


+5.498 (7.4%) 
-6.165(8.3%) 


+2.189 (2.9%) 
-2.205(3.0%) 


108.59 


+8.100(7.5% 
-9.022(8.3% 




+2.940(2.7%) 
-2.980(2.7%) 


300 


7.83 


+0.469 
-0.563 


6.0%) 
7.2%) 


+0.275 
-0.272 


3.5%) 
3.5%) 


12.09 


+0.750(6.2% 
-0.882(7.3% 




+0.399(3.3%) 
-0.397(3.3%) 


400 


1.38 


+0.070 
-0.089 


5.1%) 
6.5%) 


+0.053 
-0.052 


3.9%) 
3.8%) 


2.25 


+0.119(5.3% 
-0.148(6.6% 




+0.082(3.7%) 
-0.081(3.6%) 


500 


0.32 


+0.015 (4.6%) 
-0.019(6.0%) 


+0.014(4.3%) 
-0.013(4.1%) 


0.56 


+0.026(4.6% 
-0.034(6.1% 




+0.022(4.0%) 
-0.022(3.9%) 


niT (GcV) 


o-tot [fb] 


Scale 


PDF 


o-tot [fb] 


Scale 


PDF 


600 


90.80 


+3.945 (4.3%) 
-5.059(5.6%) 


+4.380(4.8%) 
-4.031 (4.4%) 


166.74 


+7.290(4.4%) 
-9.480(5.7%) 


+7.240(4.3% 
-6.860(4.1% 




700 


28.60 


+ 1.191 
-1.512 


4.2%) 
5.3%) 


+ 1.613 
-1.406 


5.6%) 
4.9%) 


56.04 


+2.349(4.2%) 
-3.021(5.4%) 


+2.747(4.9% 
-2.475(4.4% 




800 


9.76 


+0.391 
-0.494 


4.0%) 
5.1%) 


+0.652 
-0.546 


6.7%) 
5.6%) 


20.50 


+0.827(4.0%) 
-1.058(5.2%) 


+ 1.153(5.6% 
-0.995(4.9% 




900 


3.52 


+0.135 
-0.172 


3.8%) 
4.9%) 


+0.280 
-0.229 


8.0%) 
6.5%) 


7.97 


+0.310(3.9%) 
-0.397(5.0%) 


+0.518(6.5% 
-0.436(5.5% 




1000 


1.32 


+0.049 
-0.063 


3.7%) 
4.8%) 


+0.125 (9.5%) 
-0.101 (7.7%) 


3.24 


+0.122(3.8%) 
-0.157(4.8%) 


+0.246(7.6% 
-0.202(6.2% 




1100 


0.51 


+0.018 
-0.023 


3.6%) 
4.6%) 


+0.057(11.2%) 
-0.046(9.1%) 


1.36 


+0.049(3.6%) 
-0.064(4.7%) 


+0.120(8.9% 
-0.098(7.2% 




1200 


0.20 


+0.007 (3. 57o) 
-0.009(4.6%) 


+0.026(13.1%) 
-0.021 (10.7%) 


0.58 


+0.021(3.5%) 
-0.027(4.6%) 


+0. 060(10. 37o) 
-0.049(8.3%) 



Tabic 4: Total cross sections, at NLO+NNLL level, for the production of a heavy quark at the LHC (\/S = 7 and 8 TeV), 
including the uncertainties from scale variations and PDFs (using the MSTW2008nnlo68cl set). 



Reference [20[ uses a fixed-order approximate NNLO approach to the total inclusive cross-section, includ- 
ing the C^^'"-* terms that implement the exact 0{a'^) scale dependence. The overall uncertainty is estimated 
by just varying the scale in this framework, without accounting for the uncertainty of the finite C^^'"-* pieces, 
leading, as we argued above, to a much reduced and in our view optimistic systematics. This calculation 
has been implemented in the program H ATHOR [lOl • 

Reference [ssl ] pursues a resummation approach that shares many similarities with our work. Its authors 
resum directly the total inclusive cross-section by implementing the same anomalous dimensions and 2-loop 
Coulomb terms used here, and do not impose exact 0{a'l) scale dependence. The resummation method 



instead differs. In Ref. [55[ the so called momentum space approach of Ref. [73[ is used, which is an x- 



spacc a ppr oach, while we use an iV-space resummation, followed by a Mellin inversion. In the approach 
of Ref. [731 ■ the x-space perturbative expansion of the resummed cross section is convergent, while in our 
approach the perturbative expansion of our TV-space result is convergent, and its Mellin inversion to x s pac e 
is asymptotic. This feature has been criticized as a drawback of the Mellin space approacl^ in Ref. \m. 
However, we remind the reader that the ambiguity associated with the asymptotic nature of the Mellin 
inversion is very weak, corresponding to an effect that is suppressed more strongly than any inverse power 
of the process scale, and that in practice has totally negligible effects. Furthermore, factorization in N 
guarantees naturally momentum conservation. Although momentum conservation can be abandoned in the 
soft approximation, it was shown in [42| that it can lead to large subleading effects. We thus believe that 
the Mellin space approach is worth pursuing for this positive featurcj^. 

At the approximate NNLO order, our NNLO^ rates and scale systematics for (7^^'°^ = (Tables [1] and [5]) 
agree precisely with the equivalent results, labeled NNLOapp, in Tables 8 and 10 of [H^. After resummation, 
the differences with our work in the final predictions and theoretical systematics must be attributed to the 
different formalism. One such source of difference, for example, is that in the x-spacc resummation approach 



® In ref. FtSI it is claimed that integration over the Landau pole also arise in the computation of the Af-space resummation 
formula. We remark, however, that this pole is irrelevant for the derivation of the N space formula, that in fact has a convergent 
perturbative expansion. 

critical comparison of the x- and A'^-space method at the analytic level has been presented in ref. (74| for the case of 
Drell-Yan pair production. 



14 



additional scales are present (in our case these are only and I-Ir). Wc also note the different default values 
used for the unknown two-loop constants, which is also a reflection of the different formalisms used, since in 
the A^-space approach the resummed terms vanish for = 1 by construction|f| It is perhaps surprising that 
the largest difference among central values is observed for the Tevatron, while at the LHC central values 
are very close. This could be related to the observation made in (55| . namely that the contribution of the 
qq channel is poorl y a pproximated by the threshold expansion. Due to the dominance of this channel at 
the Tevatron, Ref. [55| argues that this could also explain why the Tevatron uncertainty does not improve 
after NNLL resummation. In all cases, the numerical differences are nevertheless consistent with the overall 
uncertainties quoted both in our work and in [55| . 

An alternative approach to the total tt cross-section has been pursued in Refs. 



[23|,126|, 1271. Like Ref. [551, 

Refs. m,!!! |27| are based on the momentum space approach. There are a number of additional differences 
between our work and these papers. They resum not the total but the differential cross-section for ti 
production. Once resummation is performed at the differential level, the differential distribution is integrated 
over phase space to obtain the total inclusive rate. In such an approach the leading terms at absolute 
threshold are correctly reproduced (see Ref. [45|) but one introduces a different set of unknown power 
corrections (as also pointed out in Ref. [IBl). Various choices for the hard scales, which become available 
when calculating differential quantities, have been explored in Refs. [H, S [ill. While the central value 



of the total cross section in Ref. [271 is slightly different from ours and from the results of Ref. [55[ , it is 



reassuring that they all still fit within the quoted uncertainty bands (the minor difference in the reference 
top mass in j^^], rutop = 173.1 GeV, has no impact in this comparison). 

Overall, the non-PDF related uncertainties in the x-space resummation approaches [27, 55 1 tend to be 
smaller than ours (by between about 25% to 40%), with the exception of the Tevatron prediction of Ref. [HBj, 
which has a corresponding uncertainty about 20% lar ger than ours. 

An approach that shares similarities with [23, 2^ 27 1 has been pursued by Ref. [Isj. In that reference 
an approximate, fixed order truncation of the differential cross section is derived and integrated over phase 
space to obtain the fully inclusive cross section. For the LHC the scale variation and central values derived 
in Ref. [1^ are similar to those of Ref. [g^ (which is about 50% smaller than our benchmark result). For the 
Tevatron the central values of these two references are also rather close, while the uncertainty of Ref. [l^ is 
much smaller than that of all other groups (it is about 60% smaller than ours); see also Refs. [tH ]. 

Similarities and differences between some of the approaches above have already been addressed in 
Refs. [t^, [H^. Such significant differences can be partially understood with the help of the discussion 
in Section [3] Overall, the large differences between central values and systematics reported in the various 
papers discussed in this section appear to be another confirmation of our conclusions about the precision 
with which the total tt cross-section is presently calculated. 
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Appendix A. Properties of the finite threshold terms 

The presently unknown constants C^^f' introduced in Eq. ^ are related to the also unknown constants 
ifj j(l) appearing in the two-loop hard matching function ([T]). With a direct calculation, and presenting 



This is to avoid introducing corrections at = 1 from a formalism that is vahd at large A'^. 
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Table A. 5: Values of Aijj in pb for = fip = m = 173.3 GeV and the MSTW2008nnlo68cl PDF set. 



Collider 


Aqg^S ^gg,S. ^gg.l 


Tcvatron 


0.3452 0.0241 0.0079 


LHC7 


1.698 4.313 1.305 


LHC14 


5.338 27.14 7.967 



directly numerical values, we obtain: 

C^!^ = -489.168+ 16 

= -1334.18+ 16 i/^^^l), 
^Si^ = -643.397+ 16 i/(^|i(l), (A.l) 

which, for the case = 0, results in the following combinations of the constants C^^f^ that enter the 



color-averaged cross-section: 

Cff^ = -1136.81 , 

Cf/^ = -489.168. (A.2) 

In an analogous way we shall define as H (1) the values of the if'^'(l) constants obtained when setting 
Cl^f = 0. We note that the constants C^^f and i?^^i(l) arc defined in different normalizations (as/(47r) 
in Eq. ([T]) and as/n in Eq. (O). 

The dependence of the cross-section on the constants cljf'^ can be estimated from: 

/^{2,0)\ /^(2,Q)\ /^(2,0)\ 



For HR = ^iF ^ rn = 173.3 GeV jSTl, and with PDF set MSTW2008nnlo68cl [56|, the values of A,jj are 
provided in Table IA.5I Clearly, reasonable variation of the unknown constants results in variation of the 
predicted cross-section by a few percent, setting an intrinsic limit to the precision of any estimate in absence 
of the full knowledge of the NNLO result. 
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